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Abstract 

We study the structural and electronic properties of isolated single-wall carbon nanotubes 
(SWNTs) under hydrostatic pressure using a combination of theoretical techniques: Continuum 
elasticity models, classical molecular dynamics simulations, tight-binding electronic structure meth- 
ods, and first-principles total energy calculations within the density-functional and pseudopoten- 
tial frameworks. For pressures below a certain critical pressure Pc, the SWNTs' structure remains 
cylindrical and the Kohn-Sham energy gaps of semiconducting SWNTs have either positive or 
negative pressure coefficients depending on the value of (n,m), with a distinct "family" (of the 
same n — m) behavior. The diameter and chirality dependence of the pressure coefficients can be 
described by a simple analytical expression. At Pc, molecular-dynamics simulations predict that 
isolated SWNTs undergo a pressure-induced symmetry-breaking transformation from a cylindrical 
shape to a collapsed geometry. This transition is described by a simple elastic model as arising 
from the competition between the bond-bending and PV terms in the enthalpy. The good agree- 
ment between calculated and experimental values of Pc provides a strong support to the "collapse" 
interpretation of the experimental transitions in bundles. 

PACS numbers: 61.46. +w, 62.25.+g, 73.22.-f 
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I. INTRODUCTION 

Since their discovery []]], carbon nanotubes have been the subject of extensive investiga- 
tion. Their variable electronic and mechanical properties M make them ideal candidates for a 
variety of applications p||. Most studies of hydrostatic pressure effects in these materials have 
focused on the analysis of structural changes in bundles 



QflBflaByflQyQ. 

The interpretation of these changes (ovalization, collapse or polygonization) seems to be 
obscured by the polydispersity of nanotube diameters in bundles jl^ lisj ]. Therefore, the 
study of the mechanical response of a single, isolated single-wall carbon nanotube (SWNT) 
under hydrostatic pressure is of significant value. 

Also of great importance is the study of changes in the electronic structure of SWNTs 
under hydrostatic pressure. Changes in the gap of semiconducting SWNTs upon different 
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brmations have been considered theoretically by many authors 
but so far the case of hydrostatic pressure has not been 
addressed. As a special motivation, recent experiments j3l in sodium dodecylsulfate (SDS)- 
wrapped SWNTs suspended in D2O show a reduction in the optical gap of all measured 
SWNTs (with varying diameters and chiralities) upon increasing pressure, a remarkable 
result that calls for theoretical interpretation. 

In this work, we address the changes in structural and electronic properties of semicon- 
ducting SWNTs under hydrostatic pressure. Section 2 describes our methodology, consist- 
ing of a hierarchy of theoretical techniques, ranging from continuum elasticity models and 
classical molecular dynamics simulations to tight-binding electronic structure methods and 
first-principles total energy calculations. In Section 3, structural and electronic properties 
in the low-pressure regime (below a critical pressure Pc) are presented. In Section 4, the 
collapse transition at Pc is described. Section 5 lists our main conclusions. 



II. METHODOLOGY 



Our first-principles calculations are based on density-functional theory (DFT) |25[ within 
the local density approximation (LDA) with the Perdew- Zunger exchange-correlation func- 



tional 



Ab initio TrouUier-Martins pseudopotentials 



are used. Calculations are 



performed using the SIESTA code |j2^], which expands the Kohn-Sham wavefunctions in a 
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linear combination of atomic orbitals (LCAO). A double-zeta plus polarization (DZP) basis 
is used. The grid cutoff representing the charge density corresponds to a plane-wave energy 
cutoff of 240 Ry. Isolated nanotubes are modeled by a periodic supercell with sufficiently 
large lateral dimensions. The irreducible Brillouin Zone is sampled with a Ixfx8 Monkhorst- 
Pack grid [2^. The effects of hydrostatic pressure P in a SWNT are simulated by imposing 
a constraint of radially directed external forces with magnitude PA, where A is the surface 
area per carbon atom. The surface area is calculated by using the external radius of the 
nanotube, i.e. the average nanotube radius plus a van der Waals "exclusion distance" of 
Ar = 1.675 A (half of the graphite interlayer distance) j^. The unit cell length along the 
axis {z direction) is adjusted so that the zz component of the stress tensor (scaled by the ra- 
tio between the supercell and nanotube cross-sectional areas) matches the target pressure P. 
This scheme provides an accurate modeling of structural and electronic properties (within a 
single-particle description) of an individual carbon nanotube under hydrostatic pressure at 
T = K. We apply it to several zig-zag SWNTs in the 5.5 A to 15 A diameter range. 

First-principles calculations are more time-consuming for chiral tubes. However, struc- 
tural properties are not strongly dependent on chirality. Therefore one can combine the 
elastic constants calculated within first-principles for zig-zag tubes with an empirical tight- 
binding (TB) model to calculate the electronic properties of chiral tubes. Although less 
reliable, the TB calculations are very useful for studying chirality trends since theyr are 
computationally less demanding. We use an orthogonal basis of one p orbital per carbon 
atom, with nearest-neighbor hopping matrix element 7 = —2.89 eV for undistorted bonds 
[2^. Under hydrostatic pressure, bond components along axial and circumferential direc- 
tions are distorted according to the radial and axial elastic constants obtained from first- 
principles. Hopping matrix elements of distorted bonds are modified according to Harrison's 
inverse- square rule jsij. 

Molecular dynamics (MD) simulations allow for the treatment of finite-temperature ef- 
fects and for a more realistic description of hydrostatic pressure by adding a pressure- 
transmitting medium. However, they are limited in accuracy owing to the use of classical 
interatomic potentials. We use a relatively well-tested and reliable potential, namely the ex- 
tended Tersoff-Brenner potential js^ to model the carbon-carbon bonding of the nanotubes. 
Hydrostatic pressure is applied to the nanotube using a method similar to that of Martonak 
et a/. The nanotube is immersed in a highly diffusive medium of particles interacting 
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via a repulsive potential. The pressure during each simulation is calculated in a box 

within the full simulation cell which contains only the pressure medium by computing the 
contributions to the stress tensor from the internal and thermal energies on the particles 
within this box. This procedure allows hydrostaticity to be easily monitored and alerts one 
to any problems which may arise from e.g. crystallization or vitrification of the medium. 



III. LOW-PRESSURE REGIME 
A. Structural properties 

For low enough pressures (below Pc), SWNTs remain cylindrical in shape. Therefore, 
structural changes under pressure are simply described by radial and axial strains, and 
e^, respectively. These strains are related to P by radial and axial elastic constants, C,. and 
Cz [SJ]! defined as Cr = —P/^r and Cz = —P/ez- Fig. [H shows Cr (squares) and Cz (circles) 
calculated by first-principles as a function of d for a few zig-zag tubes: (7,0), (8,0), (10,0), 
(11,0), (14,0), (16,0) and (19,0). As expected for a highly anisotropic material, Cr and Cz 
are not equal. The same trend was obtained in first-principle calculations in bundles 3^ . 
Moreover, notice the strong diameter dependence of both elastic constants. This dependence 
can be modeled by treating the nanotube as an elastic sheet with compliance constants Sn 
and 5*12. By using the standard definitions of Young's modulus Y = l/Su and Poisson's 
ratio u = —Su/Su, and applying geometrical considerations, we arrive at the following 
expressions fSsf: 

^ AYrAr 

' Ro{2r - vRo) ^ ' 

_ AYrAr 
Ro{Ro-2ury 

where Ro is the nanotube external radius, Ro = r + Ar. Best fits give the parameters 
Y = 1010 GPa and u = 0.242 36j. The agreement between the first-principles results and 
the elastic model is remarkable, even for small-diameter tubes. 
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FIG. 1: Elastic constants Cr (red squares) and Cz (blue circles). Lines are best fits using expressions 
dU and 0. 

B. Electronic properties 

Table m shows the calculated pressure coefficients of the LDA band gap, dEg/dP, for 
several zig-zag SWNTs. The data reveal two interesting features: (i) The sign of dEg/dP 
depends on q = {n — m) mod 3, being negative for q = 2 and positive for q = 1; (ii) The 
magnitude of dEg/dP seems to increase with diameter. Also shown in Tableware the values 
of dEg/dP calculated within TB for the same tubes, so as to compare the accuracy of the 
TB model for this quantity. The overall agreement is fairly good. The simplicity of the TB 
model allows us to use it on a wide range of diameters and chiralities. The results are shown 
in Fig. EKa). Each dot corresponds to a particular nanotube, and the calculated dEg/dP are 
plotted as a function of Eg itself. A more complete picture emerges: Values of dEg/dP seem 
to follow trends according to the specific values of {n — m), a so-called "family behavior" . 
As a guide to the eye, results for nanotubes in the same (n — m) family are grouped by color. 
The g-dependent sign oscillation found for zig-zag tubes within LDA is a manifestation of 
this family behavior. 
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TABLE I: Pressure coefficients for LDA and TB gaps of a few zig-zag SWNTs. 
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(7,0) 0.63 
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FIG. 2: a) Band-gap pressure coefficient as a function of the gap for a large number of semicon- 
ducting SWNTs. Tubes are grouped into colors according to their (n — m) family. The values of 
(n — m) for each family are also shown in the figure, b) Collapse of dEg/dP values to a single line 
when plotted against (— l)'?(icos(30), using the same color scheme as in Fig|2{a). 

Similar types of family behaviors have been found, for instance, in the study of band-gap 



changes in SWNTs under uniaxial stress 



20|,|23|. In fact, the two situations (uniaxial 



stress and hydrostatic pressure) are conceptually identical, differing only in the signs and 
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magnitudes of the resulting radial and axial strains. Therefore, one can readily adapt the 



relation for the gap shift from Gartstein et al. 



231 to the case of hydrostatic pressure: 



5Eg = + 3|^|(e^ _ e,)(-l)^cos(3e), (3) 

where ac-c is the C-C bond length and 9 is the chiral angle. The pressure coefficient then 
becomes: 

dEg _ 4|7|ac-_c oi. i f ^ ^ 



One can obtain a simpler and more useful (albeit approximate) expression by taking the 
large-diameter limit in the expressions for Cr and Cz- In this limit, both elastic constants 
decay as 1/d, and the resulting pressure coefficient is: 

dEg _ |7| 



dP 2YAr 



ac-ci2 - z/) - ^(1 + u){-iydcos{39) 



(5) 



We test this analytical expression against our numerical results by plotting dEg/dP as a 
function of {—iydcos{39) in FigEl^b). The data collapse into a straight line is excellent. 

Equation has a very clear physical meaning. Hydrostatic pressure causes an overall 
shortening of C-C bonds, therefore increasing, in a TB picture, the magnitude of hopping 
matrix elements. In a graphene sheet, this would lead to an increase in the Fermi veloc- 
ity. Because, in a simplified picture, energy gaps of semiconducting SWNTs are obtained 
by slicing the graphene bands, this would lead to an overall tendency of Eg to increase 
upon applying pressure, for all SWNTs. This is the meaning of the first term in Eq.©, 
a chirality- independent positive constant. In addition to that, hydrostatic pressure breaks 
the triangular lattice symmetry of the parent graphene sheet due to the difference in radial 
and axial strains. This leads to a relative shift of the slicing planes with respect to the 
graphene Fermi point. This shift depends on chirality: For g = 1 (2), the p 
from (closer to) the Fermi point, therefore increasing (decreasing) the gap 
meaning of the second, chirality-dependent term in Eq. (0). 

Experimentally, however, optical gaps of SDS-wrapped SWNTs in aqueous solution show 
a„ overaU *o.e.,. „.th p.e.u.e, with — d« „t ,E./,P atoost ten t.r.e. highe. tha. 
our calculated values for similar diameters [2J]. A family behavior different from the one 
described in Eq.(jSl) is found. Although excitonic effects (not included here) are crucial for 
a quantitative description of such optical experiments [22|, trends are often well described 
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by a single-particle picture. Therefore, this qualitative disagreement between theory and ex- 
periment is puzzling and it should motivate further work. It is possible that the interaction 
between the SDS micelles and the SWNTs plays an important role. Very recent experi- 
ments js^ in which strain is applied to SWNTs by both differential thermal contraction 
upon freezing the D2O solution and by hydration of a wrapping polymer provide good qual- 
itative agreement (positive and negative shifts depending on chirality) with our theoretical 
predictions. 



IV. THE COLLAPSE TRANSITION 



Several high-pressure experiments in bundles indicate some sort of symmetry-breaking 
phase transition at critical pressures ranging from 1.5 to 2.1 GPa for laser-grown tubes 
(c/ ~ 12 - 14 A) 0,y,|ll| and at 6.6 GPa for HiPCo tubes (c/ ~ 8 A) ij- One may ask the 
question: Are those symmetry-breaking transitions induced by bundling or are they related 
to an intrinsic property of the nanotubes? To investigate how isolated SWNTs respond to 
pressure, MD simulations are performed as described in Section 2. Figs. Efa) andEfb) show 
snapshots of room-temperature simulations at 3 GPa and 4 GPa, respectively, for of a chiral 
(8,7) SWNT {d = 10.3 A). At 3 GPa, the tube is cylindrical, and at 4 GPa it is collapsed 
into a fiat shape. Therefore, isolated nanotubes indeed suffer a symmetry-breaking transition 
under hydrostatic pressures. Similar conclusions have been reached by molecular mechanics 
I39I 1 and elastic models j4o|. In fact, Zang et al. 4^ suggested that a continuous change of 
cross-section shape (from circular to oval to "peanut" or collapsed) occurs upon increasing 
pressure beyond Pc- We do not observe any stable oval shape in our MD simulations. 

We also formulate a continuum elasticity model to investigate the diameter dependence 
of the critical pressure Pc- The total enthalpy per unit length of the nanotube is written as: 

"-(C-Co)^+ Ids4. + PA, (6) 



2Co' Jc 

where C is the perimeter around the tube, Cq is its value at P = 0, i? is the local curvature 
radius and A is the nanotube cross-sectional area (including the contribution from the van 
der Waals radius Ar). The first and second terms are continuum versions of the bond- 
stretching and bond-bending energies, respectively. The parameter a is estimated from the 
elasticity of graphite: a = Cwcjl = 3350 GPa. A, where c/2 is the interplanar distance of 
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FIG. 3: Snapshots of room-temperature MD simulations of a (8,7) SWNT at (a) 3 GPa and (b) 
4GPa. Large orange balls and small yellow balls represent carbon atoms in the nanotube and 
repulsive Lennard-Jones medium particles, respectively. 



graphite. The parameter (3 is extracted from first-principles calculations of strain energy of 



nanotubes with varying diameters 



41j . yielding /3 = 133 GPa.A^. The onset of instability 
of the cylindrical shape is searched by monitoring the change in enthalpy 6h when a small 
circular-to-elliptical deformation is applied to the cross-section, at different pressures. For 
P > Pc, Sh should be negative. This procedure is convenient since it allows the use of the 
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FIG. 4: Critical pressure for collapsing as a function of nanotube diameter. Open blue dots are 
the elastic-model calculated results, the blue dashed line is a guide to the eye. Full red symbols are 
experimental results. Square: Ref. up-triangle: Ref. down-triangle: Ref. ^, diamond: 

Ref. In the inset, Pc is plotted against The straight line represents the limiting 

behavior. 

well-known formulae for the ellipse perimeter, area and curvature radius. In accordance 
with Ref. we find that instability is driven by the competition between the bond- 

bending energy (that tends to keep the circular cross-section) and the PV term (that tends 
to decrease the nanotube volume by collapsing it). 

Figure IH shows the calculated values of Pc as a function of nanotube diameter. As one 
can see, P, decreases with increasing diameter, as suggested in recent theoretical works 
jisl 39 . |4q| . In the inset of Fig. HI we confirm that Pc decays as l/rf^ in the limit of 



large tubes, in agreement with the classic study by Levy on the theory of elastic rings j42 1 . 
Also shown in the figure are the available experimental data. The agreement is excellent, 
indicating that the instability of individual tubes under pressure is most likely the driving 
force for the observed phase transitions in bundles. 
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V. CONCLUSIONS 



We have investigated the structural and electronic properties of semiconducting SWNTs 
under hydrostatic pressure. For low enough pressures, nanotubes remain cylindrical and are 
described by radial and axial strains. Radial strains are larger than axial strains, a result 
that is well described by the elastic properties of a graphene sheet. Pressure coefficients 
for the single-particle band gap can be positive or negative, with strong family behavior, in 
disagreement with recent optical measurements in micelle- wrapped tubes in aqueous solution 



2J|. At a certain critical pressure Pc-, isolated SWNTs collapse from cylindrical to a flat 



shape. The good agreement between calculated and experimental critical pressures indicate 
that such a pressure- induced instability of isolated SWNTs is the driving force for the 
observed transitions in bundles. 
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